固体药柱的大变形分析

采用物质描述法，选择初始构形作为参考构形，导出以Ｋｉｒｃｈｈｏｆｆ应力和Ｇｒｅｅｎ应变表示的粘弹性大变形本构关系和有限元方程，利用ＴｏｔａｌＬａｇｒａｎｇｉａｎ增量法进行求解，取得了满意的结果。     

固体推进剂非稳态燃烧模型数值特性

在修改的KTSS燃烧模型的基础上,研究了热波燃烧模型的数值特性。在大幅度改变压力变化速率参数（β）和热释放参数（H）下,计算结果显示：1）在大β值时,燃速并未发散溢出,直到β=∞,燃速仍然是一个有限值,在一个峰值波动之后,趋向于稳态燃速。2）在增大H值时,燃速的不稳定是逐渐增加的。当H＜1时,燃速振荡后最终回到稳态值；当1＜H＜1．045时,燃速是一系列重复、有限的尖峰；但当H≥1．045时,燃速发散。     

布雷顿和卡诺热机循环性能比较(Ⅱ)定常态流变温热源循环

研究定常态变温热源热机循环性能，导出内可逆卡诺热机和布雷顿热机的最佳功率、效率关系和最大功率及相应的效率界限，并对这两种热机循环的最优性能进行了比较。理论分析表明，在相同的边界条件和热效率下，布雷顿循环的功率可以高于卡诺循环的功率，极限情况下前者是后者的两倍。对于变温热源条件，布雷顿循环主要受益于其工质与热源间有较佳的匹配。所得结果对热机工作参数和工质的最优选择有一定指导意义。     

LRE模型燃烧室燃烧过程数值模拟

建立了任意斜交曲线坐标系下液体火箭发动机（ＬＲＥ）内部工作过程的气液两相湍流化学反应流模型。应用该模型对气氢液氧模型燃烧室内部燃烧过程进行了数值模拟，得到了全流场的速度分布图形、燃气组分等值线、流场等温线、流场等压线、流场等马赫数线、壁面温度与径向平均温度曲线和燃烧效率曲线。计算结果表明：用数值模拟方法分析ＬＲＥ内部工作过程是可行的。     

一种模拟喷管二维流场的新方法

采用ＳＩＭＰＬＥ方法对固体火箭发动机喷管进行了纯气相流场的数值模拟。推导了建立在同位网格之上利用协变物理速度分量为计算变量的离散方程，采用协变物理速度分量推导压力修正方程，对密度采用一阶迎风格式插值，将ＳＩＭＰＬＥ方法扩大到计算可压流动。计算结果表明：本方法计算精度高，收敛速度快，程序编写简单。     

喷管流场和羽流场的N－S数值模拟

利用时间推进的方法，通过求解可压缩的、雷诺平均的Ｎ－Ｓ方程，并利用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数滴流模型和Ｐｒａｎｄｔｌ混合长度模型，数值模拟了飞机和火箭发动机喷管的内流和羽流，获得了喷管内流场和羽流场的参数分布，计算结果与试验结果及国外研究成果一致。网格划分采和了自适应网格技术。     

瓦状特征型面塞式喷管三维计算与实验

为了进一步了解瓦状塞式喷管的性能,采用NND差分格式求解三维N S方程和空气冷流对6单元瓦状特征型面塞式喷管进行了数值模拟和实验研究。研究模型的内喷管面积比为4,总面积比为40,设计压强比为1047。计算得到了流场马赫数和塞锥表面压强分布、喷管推力系数效率,以及不同压强比下中心平面、过渡平面和边缘平面的塞锥表面压强变化规律。计算结果与实验数据吻合得较好,效率数值最大相差1%。实验塞式喷管最大的推力系数效率为0 995,同钟型喷管相比,具有很好的高度补偿能力:从地面到高空,效率在0 93～0 995之间变化。和以前简化型面的4单元瓦状塞式喷管相比,实验和数值模拟均说明塞锥特征型面的优化设计提高了喷管性能,更充分体现了塞式喷管的高度补偿特性,可以成为未来工程应用的选择方案。     

环型超声速空气引射器零二次流流场数值研究

采用二维轴对称雷诺平均方程和标准κ-ε双方程湍流模型，数值研究了环型超声速引散器零二次流的流场结构及盲腔压强的变化，空间上采用二阶迎风格式进行耦合求解，时间上采用显示的Runge-Kutta方法进行迭代推进。结果表明，引射器几何参数不变的情况下，启动后的盲腔压强与引射气流总压之比为一常数；喷管马赫数不变情况下，喷管出口面积与混合室入口面积比越小，盲腔压强越低，扩压器性能越好，启动要求的总压越低，对超声速空气引射器的设计具有指导意义。     

二维过渡区喷流干扰流场的DSMC/EPS仿真

应用DSMC／EPS混合算法，数值模拟过渡区有RCS横向喷流干扰的二维有限平板高超声速流场，分析了流场结构、分离形态及物面气动系数分布等特性。结果表明，不同：Knudsen数下，喷流干扰流场存在差异，过渡区中的稀薄气体效应得以体现，证明了DSMC／EPS混合算法是研究稀薄，连续混合区域高超声速流动的有效方法。     

结构参数对圆转矩形喷管换热的影响

为了解结构参数对圆转矩形内喷管再生冷却换热的影响，设计了多个圆转矩形喷管，考虑了三种结构参数：转方位置、出口高宽比和出口圆角大小的影响。采用有限体积法求解三维可压缩的N-S方程对其内部流动和换热进行了数值模拟。湍流模型采用标准的k-ε双方程模型，壁面附近的流动和传热采用壁面函数法处理，速度与压力的耦合采用SIMPLE算法求解。结果表明：在型面一阶导数连续的情况下，转方位置对圆转矩形内喷管的换热影响不大；出口高宽比对圆转矩形内喷管的换热影响较大，出口高宽比不能太小，否则影响内喷管流场和换热；出口圆角大小影响内喷管周向上的温度分布，圆角太小造成周向温度分布不均匀。